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ABSTRACT 

Accretion disks with masses ~ 10~ 3 — O.IMq are believed to form during the merger 
of a neutron star (NS) with another NS and the merger of a NS with a black hole 
(BH). Soon after their formation, such hyper-accreting disks cool efficiently by neu- 
trino emission and their composition is driven neutron-rich by pair captures under 
degenerate conditions. However, as the disk viscously spreads and its temperature 
drops, neutrino cooling is no longer able to offset viscous heating and the disk be- 
comes advective. Analytic arguments and numerical simulations suggest that once this 
occurs, powerful winds likely drive away most of the disk's remaining mass. We cal- 
culate the thermal evolution and nuclear composition of viscously spreading accretion 
disks formed from compact object mergers using one-dimensional height-integrated 
simulations. We show that freeze-out from weak equilibrium necessarily accompanies 
the disk's late-time transition to an advective state. As a result, hyper-accreting disks 
generically freeze out neutron-rich (with electron fraction Y e ~ 0.2 — 0.4), and their 
late-time outflows robustly synthesize rare neutron-rich isotopes. Using the measured 
abundances of these isotopes in our solar system, we constrain the compact object 
merger rate in the Milky Way to be < IO^A/^o/O-IMq) -1 yr _1 , where M dj0 is the 
average initial mass of the accretion disk. Thus, either the NS-NS merger rate is at 
the low end of current estimates or the average disk mass produced during a typical 
merger is <C O.IMq. Based on the results of current general relativistic merger simu- 
lations, the latter constraint suggests that prompt collapse to a BH is a more common 
outcome of NS-NS mergers than the formation of a transient hyper-massive NS. We 
also show that if most short duration gamma-ray bursts (GRBs) are produced by 
compact object mergers, their beaming fraction must exceed /& ps O.^M^o/O.IMq), 
corresponding to a jet half-opening angle ;> SO^M^o/O.IMq) 1 / 2 . This is consistent 
with other evidence that short duration GRB outflows are less collimated than those 
produced in long duration GRBs. 

Key words: nuclear reactions, nucleosynthesis, abundances accretion disks — black 
hole physics — gamma rays: bursts — neutrinos 



1 INTRODUCTION 

Massive, compact accretion disks are thought to form in 
a number of astrophysical events, including the merger of a 
black hole and a neutron star (BH-NS; Lattimer & Schramm 
1974, 1976), the merger of a double NS binary (NS-NS; Eich- 
ler et al. 1989), the collapse of a rapidly- rotating stellar 
core (Woosley 1993), the accretion- induced collapse (AIC) 
of a white dwarf (WD) to a NS (Dessart et al. 2006, 2007), 
and the AIC of a NS to a BH (Vietri & Stella 1998, 1999). 



* E-mail: bmetzger@astro.berkeley.edu 



These disks are termed "hyper-accreting" due to their large 
accretion rates of up to several solar masses per second. 
Rapid accretion onto a BH following the collapse of a mas- 
sive star is a popular model for the central engine of long- 
duration gamma-ray bursts (GRBs) (MacFadyen & Woosley 
1999) due to their association with star forming regions and 
Type Ibc supernovae (SNe) (Woosley & Bloom 2006). Short- 
duration GRBs, on the other hand, originate from a more 
evolved progenitor population (e.g., Berger et al. 2005) and 
may instead result from compact object (CO) mergers (e.g., 
Lee & Ramirez-Ruiz 2007), or the AIC of a WD (Metzger 
et al. 2008a) or a NS (MacFadyen et al. 2005). 

Hyper-accreting disks cool via neutrino emission when 
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the accretion rate is sufficiently high (M J> 0.01 — O.IMq 
s _1 ; Popham et al. 1999). This allows the disk midplane to 
become dense and electron degenerate (Chen & Beloborodov 
2007). Because positron captures (e + + n — > p + v E ) are sup- 
pressed relative to electron captures (e~ + p — > n + f e ) un- 
der degenerate conditions, neutrino-cooled disks are driven 
neutron-rich (Pruet et al. 2003; Beloborodov 2003), i.e., to 
an electron fraction Y e = n p /(n p +n n ) < 0.5, where n p and 
n n are the proton and neutron density, respectively. 

If a portion of this neutron-rich material escapes the 
disk, several observable consequences may result (e.g., De- 
rishev et al. 1999). In particular, although highly relativis- 
tic outflows are required to produce GRBs, slower outflows 
from larger radii in the disk are probably at least as common 
(e.g., Pruet et al. 2004). These dense outflows may synthe- 
size heavy, neutron-rich isotopes as they expand away from 
the midplane and cool. Such outflows are difficult to detect 
directly because neutron-rich isotopes have short half-lives 
(~ seconds; Freiburghaus et al. 1999) and are thus unlikely 
to power a bright SN-like transient via their radioactive de- 
cay (although a dimmer, shorter-lived transient may be pro- 
duced; Li & Paczynski 1998; Kulkarni 2005; Metzger, Piro, 
& Quataert 2008, hereafter MPQ08). 

Neutron-rich outflows can also be probed indirectly 
via their effect on the chemical evolution of the Galaxy. 
For instance, moderately neutron-rich outflows with Y e « 
0.35 — 0.4 produce isotopes (e.g., 79 Br and 78 Se) that are 
rare in our solar system (Hartmann et al. 1985; Woosley 
& Hoffman 1992). Hence, the measured abundances of these 
isotopes place stringent constraints on the rate at which low- 
Ye material is ejected into the interstellar medium (ISM) . In 
turn, if the amount of neutron-rich material ejected in a 
given event can be estimated, the true rate and beaming 
fraction (or jet opening angle) could in principle be con- 
strained in models that associate these events with GRBs 
(e.g., Woosley & Baron 1992; Fryer et al. 1999). Constraints 
on the rate of CO mergers are particularly interesting be- 
cause they are the primary target for km-scale gravitational 
wave detectors such as LIGO and VIRGO, and rates derived 
from known merging systems (Kalogera et al. 2004; Kim et 
al. 2006) and population synthesis models (e.g., Belczynski 
et al. 2006) remain highly uncertain. 



1.1 Summary of Previous Work 

Most previous efforts to determine the composition of out- 
flows from hyper-accreting disks have focused on winds 
driven by neutrino irradiation from the inner disk (Levin- 
son 2006; Surman et al. 2006, 2008; Barzilay & Levinson 
2007; Metzger et al. 2008b), in analogy to those from proto- 
neutron stars following successful core-collapse SNe (e.g., 
Burrows, Hayes, & Fryxell 1995). However, the electron frac- 
tion at the base of neutrino-driven outflows is typically not 
preserved as matter accelerates to large radii because neu- 
trino absorptions (p + v e — > n + e + and n + u e — > p + e~) 
drive Y e to the equilibrium value 

y: ^[l + iL.A^/L.A^W 1 (i) 

given by the properties of the neutrino radiation field, where 
L Vl ,ILv P and (e Vc ) / '(ep e ) are the mean v^jv^ luminosities and 
energies, respectively, from the central NS or accretion disk 



(Qian et al. 1993). When the accretion rate is very high 
(M > Mq s _1 ) and the inner disk is optically thick to neu- 
trinos, Y e " can itself be <C 0.5 because the z> e 's originate from 
regions of higher temperature than the zVs (e.g., Surman et 
al. 2006). However, a disk with a fixed initial mass (such as 
is produced by a CO merger) spends very little time (if any) 
accreting at such a high rate; most neutrino-driven mass loss 
from a viscously-spreading disk occurs at later times, when 
the entire disk is optically thin and Y" ~ 0.5 (see Fig. 9 of 
MPQ08). Thus, while early-time neutrino-driven outflows 
may produce modest amounts of r-process ejecta, their to- 
tal yield in the context of CO mergers is probably insuffi- 
cient to contribute appreciably to the Galactic abundance 
of neutron-rich isotopes (e.g., Surman et al. 2008). However, 
as we now discuss, later stages in the disk's evolution, when 
neutrino irradiation is relatively unimportant, are likely to 
eject even larger quantities of neutron-rich material. 

In a recent paper (MPQ08), we used a one-zone (or 
"ring") model to study the evolution of hyper-accreting 
disks. Guided by analytic work (Blandford & Begelman 
1999) and numerical simulations (e.g., Stone & Pringle 2001; 
Hawley, Balbus, & Stone 2001; Hawley & Balbus 2002) 
which show that radiatively-inefncient accretion drives pow- 
erful outflows, we argued that the majority of the disk be- 
comes unbound soon following the disk's transition to an ad- 
vective state at late times. Since ejection occurs on approx- 
imately the dynamical timescale (which is faster than weak 
interactions at this stage), these "viscously-driven" outflows 
maintain the electron fraction of the disk midplane, unlike 
the neutrino-driven outflows at earlier times. Thus, in order 
to determine what isotopes are ultimately synthesized, the 
electron fraction in the disk must be known at late times 
when the disk becomes advective. 

1.2 Outline of this Paper 

In this paper we determine the composition of late-time out- 
flows from hyper-accreting disks by calculating the disk's 
midplane composition as the disk viscously spreads and falls 
out of weak equilibrium. We concentrate in particular on BH 
accretion following a CO mergeiQ (either NS-NS or NS-BH 
models), where a fixed initial disk mass is a reasonable ap- 
proximation. In 3D we use a one-zone model similar to that 
in MPQ08 to calculate Y e (t) for a wide variety of plausible 
initial disk masses and sizes. Although the one-zone model 
in SJ21 captures the basic evolution of hyper-accreting disks, 
the subsequent nucleosynthesis depends sensitively on the 
freeze-out electron fraction Y/, and so a more detailed calcu- 
lation is warranted. Furthermore, a one-zone model cannot, 
by construction, address the possibility that different an- 
nuli in the disk may freeze out with different values of Y/ . 
Therefore in §3] we present one-dimensional (ID), height- 
integrated calculations of the evolution of CO merger disks 
and their composition, and use them to determine the elec- 
tron fraction mass distribution M(Y/) at late times. We find 
that under most conditions the majority of the mass freezes 

1 Although we center our discussion on accretion following CO 
mergers, most of our conclusions would apply equally to the AIC 
of a NS to a BH, which may produce a disk with similar initial 
properties (Shibata 2003; Shapiro 2004). 
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out neutron-rich with YJ ~ 0.2 — 0.4. In Sj4]we summarize 
our results and use them to constrain the CO merger rate 
and the beaming fraction of short GRBs. 



2 ONE ZONE MODEL 

In this section we present calculations of the electron frac- 
tion at freeze-out YJ using a one-zone "ring" model similar 
to that presented in MPQ08. The simplicity of this model 
allows us to efficiently explore a wide parameter space of 
initial disk mass Md,o and radius rd,o, and it also provides 
a useful point of comparison for our ID calculations in Sj3] 
The initial conditions and relevant equations are presented 
in £|2 . 1 1 In £12.21 we present our results and show, using a 
simple analytic argument, that a moderately neutron-rich 
freeze-out (i.e., Yj < 0.5) is generically expected, relatively 
independent of the details of how the disk viscously spreads. 



2.1 Equations and Initial Conditions 

Since the accretion disks produced during CO mergers are 
created from tidally disrupted NS material, we take the ini- 
tial electron fraction to be Y e ° = 0.1, which is characteris- 
tic of the inner neutron star crust (e.g., Haensel & Zdunik 
1990a, b; Pethick & Ravenhall 1995). In most of our calcula- 
tions we take the BH mass to be Mbh = 3Mq . The mass and 
radial profile of the remnant disks formed from CO merg- 
ers are uncertain theoretically because they depend on the 
unknown supranuclear-density equation of state and general 
relativistic effects, which are now being explored in merger 
simulations (e.g., Shibata & Taniguchi 2006). Disk masses 
M d ~ 10 -3 -0.3M Q with characteristic sizes ~ 10 6 -3 x 10 7 
cm appear typical (e.g., Rasio et al. 2005; Oechslin & Janka 
2006; Shibata & Taniguchi 2006). 

The equations and assumptions employed in our one- 
zone model closely follow those in MPQ08. To provide a 
brief summary, our calculation follows the thermal and vis- 
cous evolution of the disk radius rd(t) (the "ring") that con- 
tains the majority of the mass, controls the accretion rate 
M onto the central object, and moves outwards as the disk 
accretes in order to conserve the total angular momentum 

1/2 

J oc Md,r d . Our energy equation includes viscous heat- 
ing and the dominant neutrino opacities. The model is cali- 
brated so as to reproduce the exact 5-function solution to the 
viscous diffusion equation (eq. [10]) at late times. Although 
we refer the reader to MPQ08 for details, we include here a 
brief discussion of changes and additions we have made and 
highlight the equations most important for this work. 
The electron fraction is evolved using 



— (A e + n + *u e n) 



Y e 



(^) 



(2) 



where d/dt is a Lagrangian derivative, A e - p /A e + n are the 
pair capture rates (see Beloborodov 2003, eqs. [6], [7]), X{ 
is the free nucleon mass fraction in nuclear statistical equi- 
librium (NSE) (from, e.g., Woosley & Baron 1992), and 



are the neutrino capture rates (e.g., Qian & Woosley 1996) 
due to irradiation from the inner accretion disk (see the dis- 
cussion below). In equation ©, vN stands for either v e n 
or v c p, Xm is the corresponding proton or neutron mass 
fraction, L„ = L52IO 52 ergs s _1 is the neutrino luminos- 
ity, (e„) = 10(eio) MeV is the mean neutrino energy, and 
cr„jv — 9 x 10 _44 (e 2 ) MeV _2 cm 2 is the neutrino capture 
cross section. The factor (1 + t v m) accounts for the possibil- 
ity that the absorbing annulus may become optically thick, 
where t v n = cj v nXn prd/2m p is the neutrino optical depth 
and p is the density; our calculation is fairly insensitive to 
this prescription for when t v > 1, however, because t v is 
generally < 1 at freeze-out. We neglect the proton- neutron 
mass difference A » 1.3 MeV and the effects of electron de- 
generacy in calculating o v n because the energies of the neu- 
trinos that dominate the heating ((e„) ~ 10 MeV) greatly 
exceed A and the electron Fermi energy, respectively. We 
assume relativistic e~/e + pairs when calculating the elec- 
tron chemical potential; we have verified that this is a good 
approximation by checking that including the effects of ar- 
bitrary e~/e + energies has no significant effect on the value 
of Yj . NSE is a generally a good assumption in calculat- 
ing Xf because the entropy in the disk is sufficiently low 
that a— particles form while the midplane temperature is 
still high (T £ 0.5 MeV). 

The dominant process that sets Y e during the thin 
disk phase is pair capture on free nuclei (e.g., Pruet et 
al. 2003). However, neutrino absorptions become important 
near freeze-out and must be included in a detailed calcula- 
tion (see Fig. 0. The neutrino absorption rate A„jv is dom- 
inated by the high energy neutrinos that are radiated from 
the smallest radii in the disk. We calculate L v and (e„) by 
extending our solutions down to the innermost stable or- 
bit r*, assuming that M is constant with radius interior to 
Td- The neutrino flux is reasonably approximated as a "light 
bulb" at the origin in equation (J3J) because the neutrino cap- 
ture rates only become comparable to the pair capture rates 
when the disk has spread to a radius that is much larger 
than that of the central BH. 

Since neutrino captures only affect Y e appreciably at 
times when the inner disk is optically thin to neutrinos, we 
assume that L„ c ~ L^ e and (e Ve ) ~ (ep e ) (see the discussion 
in §5.3 of MPQ08 for why this is a good approximation). 
Thus, the effect of neutrino captures is to drive Y e towards 
Y^? ~ 0.5 (see eq. Q~]). The effect of a-particle formation 
(i.e., Xf — > in eq. [2]) would also be to drive Y e — > 0.5 (the 
"alpha-effect"; e.g., Fuller & Meyer 1995); however, this ef- 
fect is unimportant in our calculations because weak freeze- 
out generally precedes a-particle formation. 

One difference relative to MPQ08 is that here we take 
the viscous stress to be proportional to just the ion gas pres- 
sure -Pgas, as opposed to the total pressure Ptot, which also 
includes radiation pressure from photons and e~/e + pairs. 
In other words, we take the kinematic viscosity to be 



v = aP gas /pf2, 



(4) 



4vrr 2 (£„)(! + t vN ) 



rg(l + t„ at) 



(3) 



where a is a dimensionless constant and SI is the Keplerian 
rotation rate. At early times in the evolution of the disk (in 
particular when it is neutrino-cooled), the disk is primarily 
supported by gas pressure and so using either v oc Ptot or 
v oc Pgas gives similar results. However, as the disk becomes 
thick at late times and Y e freezes out, the disk becomes 
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Figure 1. Contours of the final electron fraction Yj following 
weak freeze-out in CO merger disks as a function of the initial 
disk mass q and radius q for two values of the viscosity, 
a = 0.3 (top) and a = 0.03 (bottom). The initial electron fraction 
in all models is taken to be Yj = 0.1. 

dominated by radiation pressure, with Ptot /Pgas ~ 2 — 3 
during freeze-out (see Fig. [2}. The primary reason that we 
use Pg a s instead of P to t in equation @ is to avoid the clas- 
sic Lightman-Eardley (1974) viscous instability, which we 
otherwise find develops in our full ID calculations described 
in 331 Because Ptot /Pgas is modest (< 3) during the times 
that matter to our results, we choose a ~ 0.3 and 0.03 as 
fiducial values for our calculations, to roughly bracket the 
value of a ~ 0.1 motivated by a number of astrophysical 
observations (King et al. 2007). For typical disk parameters, 
we find that taking v oc P gas instead of v oc P tot in our one- 
zone models results in a modest (~ 20%) increase in the final 
electron fraction. This, together with our inclusion of neu- 
trino absorptions in evolving Y e , accounts for the fact that 
the final electron fractions YJ presented here are somewhat 
larger than those presented in MPQ08 (their Fig. 9). 

As an improvement over MPQ08, our calculation of the 
cooling rates due to e~/e + captures on free nuclei now in- 
cludes the effects of arbitrary electron energy and degener- 
acy. In addition, we include heating from the absorption of 
neutrinos that are radiated from the accretion disk at small 
radii (see eq. [14]). Finally, as in MPQ08 we apply a "no 
torque" boundary condition at the inner radius r* = 10 6 
cm, which corresponds to the radius of the innermost stable 
orbit around a rapidly rotating Kerr BH (a = 0.9). 

2.2 Results 

Figure [1] shows contour plots of the freeze-out electron frac- 
tion Y/ as a function of the initial disk mass Md,o an d radius 
rd,o using two different values of the viscosity, a = 0.03 and 
0.3. 

The basic features of Figure [T] can be understood as fol- 



lows. Sufficiently massive, compact disks (upper left corner 
of Fig. [TJ come into pair capture-equilibrium and obtain an 
electron fraction that is relatively independent of the initial 
value Y° = 0.1, while low mass, more extended disks (lower 
right corner) retain their initial composition. As discussed 
below, a disk enters (3— equilibrium if and only if it can cool 
efficiently via neutrino emission (see eq. [9]). By equating 
the radius where neutrino cooling balances ~ 1/2 of the vis- 
cous heating in a self-similar solution for a thick disk, we 
estimate analytically that disks that satisfy 

far »° < 5 » 

are sufficiently compact to become neutrino-cooled and 
neutron-rich at early stages in their evolution (see also 
MPQ08 eqs. [19]- [21]). Equation (J5]l does a reasonable job 
of reproducing the parameter space in Figure [T] where Y - / 
differs significantly from its initial value. For instance, equa- 
tion ([5]) shows that lower values of a allow less massive disks 
of a fixed size to enter /3-equilibrium; this is consistent with 
the wider parameter space with Yj > Y° = 0.1 at low a. 

To illustrate the process of weak freeze-out explicitly, 
Figure [2] (left panel) shows the evolution of the electron frac- 
tion Y e (t) for a disk with M dfi = O.IMq and r d , = 3 x 10 6 
cm. We also show the equilibrium electron fraction Y e e9 (ob- 
tained by setting the right hand side of eq. [2] to zero), the 
ratio of the disk scaleheight to the radius H/ra, the ratio 
of the gas pressure to the total pressure P ga s/Ptot, and the 
ratio of the electron chemical potential to the disk tempera- 
ture fj, e /kT (i.e., the degeneracy parameter). Also shown for 
comparison in the right panel of Figure are the weak in- 
teraction rates from equation ((2]), normalized to the viscous 
timescale i v j sc = a~ 1 Sl~ 1 (///r [ j) _2 . 

Figure [2] shows that at early times, when the disk is 
neutrino-cooled and H/ra ~ 1/3, pair captures are rapid 
compared to the viscous timescale. Thus, the disk enters 
equilibrium with an electron fraction Y e « Y/ 9 that 
is <C 0.5 because positron captures are suppressed rela- 
tive to electron captures under the degenerate conditions 
(fi e S> kT) in the midplane. At later times (t ~ 0.1 s) the 
disk thickens to H/r<i ~ 0.5 and becomes non-degenerate 
(fi e <C kT) and radiation pressure-dominated. At this point 
Y e cq rises to J> 0.5 because positron captures under non- 
degenerate conditions now become energetically favored over 
electron captures due to the proton-neutron mass differ- 
ence. Note, however, that as the disk thickens weak inter- 
actions become slow compared to the viscous timescale and 
Y e freezes out at the value Yj » 0.43. Although (non-local) 
neutrino absorptions arc unimportant relative to electron 
captures at early times, their rate A„jv is comparable to the 
pair capture rates near freeze-out and must be included to 
accurately determine Y[ . 

Figure [T] shows that moderately neutron-rich freeze-out 
(Yl ~ 0.3 - 0.5) is a common feature of accretion disks 
formed from CO mergers for a wide range of initial condi- 
tions. However, contrasting the disk's neutron-rich equilib- 
rium state (Y e — Ye q ~ 0.1) at early times with its final, 
non-degenerate state that actually favors a proton-rich com- 
position in equilibrium (i.e., Y e e9 > 0.5), it is a priori unclear 
why y e freezes out with a value Y/ ~ 0.3 — 0.5 between these 
two extremes. 

To understand why this is the case, note that an approx- 
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Figure 2. The process of weak freeze-out in our one-zone model of a viscously-spreading, hyper-accreting disk. The left panel shows 
the time evolution of the electron fraction Y e (solid line), disk thickness H/r^ (dashed line), degeneracy parameter fi e /kT (triple-dot- 
dashed line), and the ratio of gas to total pressure Pgas/Ptot (dot-dashed line), calculated for a = 0.3 and for a disk with initial mass 
M^ q = O.IMq and radius r^o = 3 X 10 6 cm. Also shown with a dotted line is the equilibrium electron fraction Y^ q obtained by setting 
the right hand side of equation J2J to zero. The right panel shows the rates of the weak interactions that modify Y e (t), normalized to the 
viscous timescale t v i sc . The electron, positron, and neutrino capture rates are denoted by X ep , A e + n , and A„jv> respectively. Note that 
Y e freezes out of equilibrium at the same time that the disk becomes advective (H/r^ > 0.5) and non-degenerate (/i e < kT). Since Y^ q 
rises to > 0.5 as the disk becomes non-degenerate, the final composition is only moderately neutron-rich, with a final electron fraction 

Yj ~ 0.43. 



imate condition for weak freeze-out is that the pair capture 
timescale becomes longer than the viscous time tvisc Since 
the disk becomes non-degenerate around freeze-out, electron 
and positron captures occur at a similar rate (e.g., Qian & 
Woosley 1996) 

XeN ~ 0.45Tm cV s -1 , (6) 

where T = TMovMeV and eN represents either e~p or e + n. 
The associated URCA cooling rate per nucleon is given by 

QeN = MA eJV ~ 2.3T& eV MeVs" 1 , (7) 

where (e„) w 5.04fcT is the mean neutrino energy released 
per capture. The viscous heating rate per nucleon at r 3> r* 
is 

Qvisc = ^vtfmn « ^amNr 2 fl 3 (^-^j , (8) 

where m,N is the nucleon mass. If we assume that the disk 
is supported primarily by gas pressure^ then the midplane 
temperature is kT « m-^r 2 ^ 2 (H / r) 2 . Combining this with 
equations (0 — JgJ we find that 

A e JV ^visc — 

Equation © shows that once neutrino cooling no longer off- 
sets viscous heating (i.e., Q v i sc > Qejv) and the disk thickens 
to H/r ~ 1, weak freeze-out (i.e., A e jvt v isc < 1) necessarily 
results. Physically, this occurs because the disk is cooled by 

2 As discussed in i|2,ll and shown in Figure [2] gas and radiation 
pressure are comparable at freeze-out; the assumption that gas 
pressure dominates is, however, sufficient for the purposes of a 
simple estimate. 



e~/e + captures, the very same processes that largely con- 
trol the evolution of Y e (see Fig. [2}. This conclusion is fairly 
robust because it applies at any radius and is independent 
of the value of a. Although the absorption of neutrinos from 
the central disk threatens to complicate this simple argu- 
ment, the right panel of Figure [2] shows that the effect of 
absorptions is at most comparable to that of pair captures 
prior to freeze-out. 



3 HEIGHT-INTEGRATED MODEL 

The results of ij2] suggest that neutron-rich freeze-out is 
a common feature of viscously-spreading, hyper-accreting 
disks. However, the one-zone model makes a number of sim- 
plifying assumption that may affect the final electron frac- 
tion Yj . In this section we present a ID (height-integrated) 
model of a viscously-spreading disk that more precisely de- 
termines Yj and its distribution with mass. In H3.ll we de- 
scribe the initial conditions and relevant equations. Then we 
present our results in H3.2I 

3.1 Equations and Initial Conditions 

Our ID calculation evolves the surface density E, midplane 
temperature T, and electron fraction Y e as a function of 
radius r and time t using the 2N-RK3 (6th order space, 3rd 
order time) scheme described in Brandenburg (2001). We 
use a logarithmic radial grid that extends from just outside 
the inner edge of the disk at r* = 10 6 cm, out to a radius 
that safely exceeds the outer edge of the disk at the final 
time step (typically w 10 9 cm). 

The surface density is evolved according to the diffu- 
sion equation for an axisymmetric disk in a Newtonian 1/r 
gravitational potential (e.g., Frank et al. 2002): 
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8E 
dt 



3d_ 

r dr 



1/2 



(*,Er 1/2 ) 



— (vXr 
or 



(10) 



where the viscosity v is proportional to just the gas pressure 
(see eq. [4] and surrounding discussion). The radial velocity 
v r is not evolved explicitly but instead follows from equation 
(|10|) and mass continuity: 



-3 d 



(^Er 1/2 ) 



ErV2 Qr 

We take the initial surface density of the disk to be 



(11) 



E(r,t = 0) oc 



r<i,o 



exp 



-(2 + m) 



r<i,o 



(12) 



with the constant of proportionality set by requiring that 
the total disk mass equals Md,o- Equation (|12p concentrates 
the initial disk mass oc Er 2 about the radius rd,a, with larger 
values of the parameter m resulting in a more narrowly con- 
centrated mass distribution. As discussed in H2.ll the precise 
mass distribution of disks produced from CO mergers is un- 
certain theoretically. In most of our calculations we take 
m = 5, although our results are relatively insensitive to 
m. Due to numerical issues that arise from the exponen- 
tial drop-off of E at the outer edge of the disk, we impose 
an initial density floor that is sufficiently small that matter 
at the density floor contains a fraction <J 10 -3 of the total 
disk mass at any time; we have verified that our results are 
insensitive to the level of this floor as long as it contains 
negligible mass. We have also performed calculations using 
a power-law density distribution that concentrates the disk's 
initial mass at small radii (e.g., E oc r -3 ); we find results 
that are similar to those obtained using equation (| 12f) with 
rd,o taken near the disk's inner edge. 

We calculate the scaleheight H by assuming hydrostatic 
equilibrium in the vertical direction, i.e., we take H — c s /fi, 
where c s = (Ptot/p) 1 ^ 2 is the isothermal sound speed. We 
evolve the midplane temperature using the equation for the 
specific entropy S, which is given by 



T— = q via 
dt 



(13) 



where (d/dt) = (d/dt) + v r (d/dr), q visc = (9/4>fi 2 is the 
viscous heating rate, and q„ is the neutrino cooling rate, 
as given in Di Matteo et al. (2002) but modified to include 
the effects of arbitrary Y e and electron degeneracy (e.g., Be- 
loborodov 2003). The term 

•+ _ i _ / _ \ f^- — xf 

1v — ^eP^I'e/ 



Y e - 



+ 



Y, : 



l-Xt 



(14) 



represents neutrino heating from the inner radii of the ac- 
cretion disk, where A„ e „ and \& a p are the neutrino capture 
rates given by equation © . We not not include heating due 
to q— particle formation because we are primarily interested 
in determining the final electron fraction, and a— particles 
form only after the disk has fallen out of weak equilibrium. 

The initial temperature profile is set by requiring that 
H/r = 0.3 at all radii. An initially thick disk is physically 
motivated by the fact that CO merger disks form dynami- 
cally hot; however, because the thermal time is short com- 
pared to the viscous time, the disk's evolution quickly be- 
comes independent of the initial scaleheight anyways. 
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Figure 3. Local disk mass 7rr 2 S and the absolute value of the 
mass advection rate M = 2nrv r T, as a function of radius. The 
model assumes a = 0.3, M^ q = O.IMq and r^fl = 3 X 10 6 cm 
(with m = 5 in eq. 1121 ). Solutions are shown at t = 0.01,0.1, 
and 1 s, with later times denoted by increasingly thinner lines. 
The nodes in the lower panel occur at the stagnation point that 
separates matter accreting onto the BH (M < 0) at small radii 
from the bulk of the disk that is spreading outwards (M > 0). 



Finally, we evolve the electron fraction profile Y e (r) us- 
ing equation ([2]), with d/dt = d/dt + v r d/dr. To keep the 
required timesteps reasonable, we equate V e with its equi- 
librium value y/ 9 when the weak interaction rates greatly 

exceed the local visCOUS TcLtG 1 / £yisc ■ 

As in SJU we take the 
initial electron fraction to be Y® = 0.1 at all radii. 

The surface density across the ghost zones at the inner 
boundary is set to enforce a constant mass accretion rate 
M oc T,v at the value of M in the first active zone. We in- 
terpolate the temperature and electron fraction at the inner 
boundary. Unlike in the one-zone calculations in [JU we do 
not implement a no-torque boundary condition on the inner 
edge of the disk, in part because it leads to an unphysical 
temperature profile at small radii (in particular, T — > as 
r — > r*). Our results for Yj are not sensitive to this bound- 
ary condition. We also interpolate all variables at the outer 
boundary, although this choice has no effect on our results 
because the outer boundary is not in causal contact with 
the inner flow. We have verified that our code conserves to- 
tal mass Mtot and angular momentum J by checking that 
decreases in Mtot and J with time are compensated by their 
fluxes across the inner grid cell. 



3.2 Results 

Figures SI and [S] summarize our results for a disk with 
a = 0.3, Md,o = O.IMq, and r dfi = 3 x 10 6 cm, which we 
take as our fiducial model. Figure[3]shows the local disk mass 
7rr 2 E (top) and the absolute value of the mass advection 
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r (cm) 

Figure 4. Scaleheight H/r, the ratio of neutrino cooling g coo i = 
q~J to total heating gh ca t, and the midplane temperature T for 
the same model and times presented in Figure [3] The total heat- 
ing ghcat = 9visc + Qv consists of viscous heating and neutrino 
irradiation from the accretion disk at small radii. 



rate M = 2nr'Ev r (bottom) as a function of radius r at 
times t = 0.01,0.1, and 1 s, with later times denoted by 
increasingly thinner lines; these epochs correspond to ~ 1, 
10, and 100 times the initial viscous time at r = rd,o- 

Figure [3] shows that the disk spreads outwards in time, 
reaching ~3x 10 8 cm by t ~ 1 s. After a few viscous times, 
a constant inward accretion rate is established at small radii 
that roughly matches the outward mass advection rate of the 
bulk of the disk. The inner steady-state disk almost extends 
to the stagnation point where v r — 0, which moves outwards 
with time. 

Figure |3] shows radial profiles of the disk thickness H/r 
(top), the ratio of neutrino cooling g coo i = to total heat- 
ing ghcat = fjvisc + qt (middle), and the midplane temper- 
ature (bottom). Figure |3] illustrates that at early times the 
bulk of the disk becomes efficiently cooled by neutrinos (i.e., 
^cooi ~ gheat) and geometrically thin (H/r <C 1). At later 
times, as the disk spreads and the temperature decreases, 
neutrinos are no long able to cool the majority of the mass ef- 
ficiently (i.e., g coo i <C Qheat) and H/r increases. As discussed 
in MPQ08, the outer disk becomes thick first, and radia- 
tively inefficient conditions move inwards as M decreases. 
This behavior can be seen explicitly in Figure [4] by compar- 
ing where H/r becomes large at early and late times. Note 
that the small "bump" in H/r corresponds to low-density 
material on the very outer edge of the disk which separates 
the density floor (where the thermal time is always much 
longer than the evolution timescale) from the bulk of the 
disk. This artifact of our initial conditions has no effect on 
our conclusions. Also note that at late times (t = 1 s) a 
range of radii around r ~ 40 km has q coo i > fjWat- This re- 
gion of net cooling radiates the thermal energy carried into 
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Figure 5. Electron fraction Y e (solid line) and the equilibrium 
electron fraction Y// q (dotted line) for the same model and times 
presented in Figures [3] and [4] At early times, the disk is in 
/3— equilibrium, with Y e Y// q >C 0.5 at all radii that contain 
substantial mass. At later times, Y// q rises as the disk becomes 
advective and non-degenerate, and Y e freezes out with final values 
in the range Yj ~ 0.1 — 0.5 (see Fig.0. 



the neutrino-cooled portion of the disk from the higher en- 
tropy advective disk at larger radii. 

Figure [S] shows the evolution of the electron fraction Y e 
(solid line) and its equilibrium value Y/r q (dotted line) . At 
very early times the entire disk is in weak equilibrium, with 
Y e = Y^ q at all radii that contain substantial massif] How- 
ever, as the disk spreads and thickens, the midplane becomes 
non-degenerate, which causes Y^ q to rise. As this occurs, 
weak interactions become slow compared to the timescale 
over which the disk evolves (for the same reasons discussed 
in £|2.2[1 , and Y e begins to freeze out of equilibrium as it lags 
behind the rising value of Y// q . As with radiatively ineffi- 
cient conditions, freeze-out begins at the outer edge of the 
disk and moves inwards with time. By the final time step 
the majority of the disk mass has frozen out, with values 
of Y/ that span the range Yj ~ 0.1 — 0.5. This behavior 
is directly analagous to the freeze-out of the "ring" radius 
in the one-zone model (Fig. [2}, but now occuring in each 
annulus of the disk. 

In Figure [6] we quantify the global process by which the 
disk becomes advective and falls out of weak equilibrium by 
showing the time evolution of the total disk mass M to t, the 
total mass that is advective M t ^ v , and the total mass that 
has fallen out of weak equilibrium M t f ot . We define annuli 
that are advective and have fallen out of equilibrium as those 
that satisfy 



<?cool <J 9he a t/2 

and 



< 3/£v 



(15) 



(16) 



respectively, where A ma x is the maximum of the weak inter- 

3 At radii larger than the extent of the disk, Y e remains equal to 
its initial value at all times because these regions (which comprise 
the density floor) never enter equilibrium. 
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the total mass that has become advective because it cannot cool 0.1 0.2 0.3 0.4 0.5 0.6 

efficiently M?~? (as defined by eq. I15I ), and the total mass that Y e 
has fallen out of weak equilibrium t (as defined by eq. I160 . 



action rates in equation ([2]). Although the numerical prefac- 
tors used in equations (|15p and (|16jl are somewhat arbitrary, 
our results are not sensitive to their precise values. Figure [5] 
shows that M to t decreases with time as matter accretes onto 
the BH, with M to t oc t -1 ^ 3 at late times, as expected from 
the self-similar behavior of an advective disk with H/r ~ 
constant (MPQ08). After an initial viscous time, a compa- 
rable amount of material is advective and out of equilibrium; 
M t a d t v and M/ ot remain fairly constant at ~ 2 x 10~ 2 M Q un- 
til the majority of the disk becomes advective and freezes 
out at late times. 

As discussed in fjl] once a portion of the disk becomes 
advective, a viscously-driven outflow likely drives away a 
substantial portion of its mass (e.g., Blandford & Begelman 
1999). In our calculations we have neglected the effects of 
such a mass sink on the evolution of the disk. Figure[U]shows 
that this approximation remains reasonable until t ~ 0.03 — 
0.1 s because prior to this point the majority of the disk 
remains neutrino cooled. Mass and angular momentum loss 
to a wind leads to a more rapid decline in M, which further 
speeds up the onset of the advective phase and weak freeze- 
out at smaller radii (MPQ08). 

In Figure [7] we show histograms of mass as a function of 
electron fraction M(Y e ) at the times t = 0.03 s and t = 0.3 
s. The solid and dashed lines are the matter that is ad- 
vective and out of j3— equilibrium, respectively (as defined 
by eqs. [15] and [16]). These distributions typically corre- 
spond to the same matter because annuli fall out of equi- 
librium as they become advective (see eq. [§])■ We show the 
composition at these two particular times because they cor- 
respond to epochs when the disk has become w 50% and 
90% advective, respectively (Fig. [BJ. Thus, these distribu- 
tions likely represent typical values of the electron fraction 
in the viscously-driven outflows that dominate subsequent 
mass loss from the disk. 

It is useful to compare the results in Figure [7] with the 
one zone calculation from for the same initial disk pa- 
rameters (see Fig. [!}. When the disk is 90% advective at 
t — 0.3 s, the height-integrated calculation finds a final elec- 



Figure 7. Amount of mass with a given electron fraction M(Y e ) 
for the model with Mdfi = O.IMq, r^o = 3 X 10 6 cm, and a = 
0.3. The solid and dotted lines include, respectively, mass that 
has become advective (eq. 1151 ) and that has fallen out of weak 
equilibrium (eq. 1161 ). The times t = 0.03 s and 0.3 s correspond 
to when the disk is 50% and 90% advective, respectively. 



tron fraction distribution centered around a mass-averaged 
value (Y/) pa 0.29, which is somewhat lower than the value 
of Yl = 0.43 given by the one-zone model. Indeed, in most 
of our calculations we find that Y/ calculated using the one 
zone model tends to slightly exceed (Y/) obtained from the 
full ID calculations. 

In Figures [8] and [9] we show histograms similar to Fig- 
ure [7] but for solutions calculated with a lower viscosity 
(q = 0.03) and a lower initial mass (Md,o = O.OIMq), re- 
spectively. The times chosen correspond to epochs when the 
disk is 50% and 90% advective. The mass-averaged electron 
fraction when the disk is 90% advective is (Y/) ~ 0.3 for 
both the lower a and lower Md,o models. The similarity in 
the distribution of Yj for a = 0.3 and a = 0.03 in Figures 
[7J and [5] supports our argument that the freeze-out process 
is relatively independent of the details of how the disk vis- 
cously spreads (see £ I2.2[) . The lower disk mass case (Fig. [9} 
is notable because a significant fraction of the disk's mass 
freezes out with Y e < 0.2, and such low- Ye material may pro- 
duce third-peak r-process elements (Hoffman et al. 1997). 

We have performed freeze-out calculations for a num- 
ber of other disk parameters. In Table [1] we summarize our 
results for (Y/) and the fraction / a d v of the initial disk mass 
that remains when the disk is 90% advective. Although dif- 
ferent initial mass distributions and viscosities give some- 
what different Y/ distributions, the neutron-rich freeze out 
of ~ 20 — 50% of the disk's original mass is a generic property 
of the disks created during CO mergers. 
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Table 1. Properties of Freeze Out in ID Calculations 



t = 1.3 s 
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Figure 8. The same as Figure [7] but for a = 0.03. The times 
shown, t = 1.3s and 8 s correspond to when the disk is 50% and 
90% advective, respectively. 
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t = 0.02 s 
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Figure 9. The same as Figure[7] but for M d>0 = O.OIMq. 



4 DISCUSSION 

Accretion disks formed from CO mergers that are suffi- 
ciently massive and compact to satisfy equation (J3J) become 
neutrino-cooled, degenerate, and neutron-rich early in their 
evolution. Eventually, as the disk accretes and viscously 
spreads, neutrino cooling is no longer able to offset viscous 
heating and the disk becomes advective. At this point sev- 
eral important changes occur nearly simultaneously: the disk 
becomes geometrically thick (H/r ~ 1), radiation pressure- 



M rf o (M(7)) 


Tri n (cm) 




M B h (Mq) 




f(b) 
J adv 


0.1 


3 x 10 6 


0.3 


3 


0.29 


0.20 


- 


3 x 10 6 


0.03 


3 


0.30 


0.15 




6 x 10 6 


0.3 


3 


0.26 


0.31 




1.2 x 10 7 


0.3 


3 


0.21 


0.40 




6 x 10 6 


0.3 


10 


0.37 


0.28 


0.01 


3 x 10 6 


0.3 


3 


0.28 


0.28 




3 x 10 6 


0.03 


3 


0.34 


0.21 




1.2 x 10 7 


0.3 


3 


0.21 


0.40 


0.001 


3 x 10 6 


0.3 


3 


0.28 


0.39 




1.2 x 10 7 


0.3 


3 


0.17 


0.72 



(a) Mass-averaged freeze out electron fraction when the disk is 
90% advective; (b) Fraction of the initial disk mass that remains 
when the disk is 90% advective. 

dominated, and non-degenerate at the same time that the 
equilibrium electron fraction Y/ 9 rises from <C 0.5 to > 0.5 
and Y e freezes out (Fig. [2] see also Beloborodov 2008). Weak 
freeze-out necessarily accompanies the disk's advective tran- 
sition because e - /e + captures are primarily responsible for 
both cooling the disk and setting its electron fraction (see 
eq. [9] of N2.2p . Because Y^ q is rising as freeze-out occurs, 
moderately neutron-rich freeze-out with Yj ~ 0.2 — 0.4 is 
a generic feature of accretion disks formed from CO merg- 
ers. Our calculations in iJ5]and 32U ar g ei y confirm this basic 
conclusion (see Table [TJ. 

Neutron-rich freeze-out in hyper-accreting disks is use- 
fully contrasted with the proton — rich freeze-out (Y/ ~ 
0.88) in the very early universe. Big bang nucleosynthe- 
sis (BBN) occurred on a timescale of minutes under non- 
degenerate and highly radiation-dominated conditions (en- 
tropy S ~ 10 10 fcs baryon -1 ). By contrast, hyper-accreting 
disks freeze out on a timescale ~ 0.1 — Is (depending on a) 
under moderately degenerate conditions, and with compara- 
ble gas and radiation pressure, i.e., S ~ 3 — 20 ks baryon -1 . 
Hyper-accreting disks freeze-out neutron-rich because they 
start in /3— equilibrium under degenerate conditions, while 
the early universe never had Y e < 0.5 because degeneracy 
effects were never important. 

A further important difference is in the production of 
heavy elements. Once the disk becomes advective and freezes 
out, a viscously-driven wind likely unbinds most of its re- 
maining mass; as these outflows expand away from the mid- 
plane and cool, heavy elements will be formed. Due to the 
"deuterium bottleneck," BBN produced almost no elements 
heavier than He. By contrast, the late-time outflows from 
hyper-accreting disks possess modest entropiefl and are thus 
generally in NSE when a— particles form. Even an a— rich 
frcczc-out is unlikely (Hoffman & Woosley 1992) because, 
under dense (low entropy) neutron-rich conditions, q's burn 
efficiently via the reaction sequence 4 He(a n,7) 9 Be(a,n) 12 C 



4 The entropy in the outflow is probably only a few ks baryon -1 
larger than in the disk midplane because only a fraction of the 
gravitational binding energy is required to drive an outflow (ad- 
vective disks are only marginally bound). 



© ???? RAS, MNRAS 000,ITJtT2l 



10 B. D. Metzger, A. L. Piro, E. Quataert 



(Delano & Cameron 1971). As discussed further in S|4. 1 1 
the abundances of the elements synthesized are thus ap- 
proximately determined by neutron-rich NSE (Hartmann et 
al. 1985). 

It is also important to contrast the site of nucleosynthe- 
sis that we have introduced here with outflows driven from 
the disk by neutrino heating (e.g., Levinson 2006; Barzi- 
lay & Levison 2008; Metzger et al. 2008b). Neutrino-driven 
outflows generally occur at early times or small radii, when 
the disk midplane is still in /3— equilibrium and the irradiat- 
ing neutrino flux is substantial. The outflow's final electron 
fraction in this case is set by an equilibrium with neutrino 
absorptions (Y e ~ Y£ in eq. [T]) and is not directly related 
to F e in the midplane (e.g., Surman & McLaughlin 2004). 
By contrast, the outflows considered here occur after the 
midplane itself has fallen out of weak equilibrium, and so 
the disk's electron fraction YJ is retained by the outflow. 
While neutrino absorptions in the disk are somewhat impor- 
tant near freeze-out, at no point does the midplane electron 
fraction enter equilibrium with the neutrinos, as generically 
occurs in neutrino-driven outflows (Metzger et al. 2008b). 
Instead, the outwardly-spreading disk provides a protective 
"cocoon" that transports low- Ye material (kept neutron-rich 
by degenerate electron capture) away from the deneutron- 
izing flux of neutrinos from the central disk. Furthermore, 
the viscously-driven outflows that we have focused on carry 
away a substantial fraction of the disk's initial mass; because 
neutrino-driven outflows likely eject less mass (MPQ08), the 
total nucleosynthetic yield from CO mergers is probably 
dominated by the viscously-driven outflows at late times. 

Finally, although we believe that the disk is largely un- 
bound by viscous heating once it becomes advective (Bland- 
ford & Begelman 1999; Stone & Pringle 2001), if nothing else 
the energy released by a— particle (or even heavier element) 
formation is sufficient to unbind the disk at late times (Lee 
& Ramirez-Ruiz 2007; MPQ08). Since a— particle formation 
occurs soon after the disk becomes advective, our conclu- 
sions would be qualitatively unchanged if nuclear (rather 
than viscous) energy powers the late-time outflow. In this 
case the formation of heavy elements would begin in the disk 
midplane, but the final nucleosynthetic yield would probably 
be similar. 

4.1 Implications 

As we have argued, moderate entropy, neutron-rich outflows 
appear to be a common property of accretion disks formed 
from CO mergers at late times. Although detailed calcula- 
tions of the dynamics and nucleosynthesis in these outflows 
need to be performed, the late-time conditions in the disk 
favor the production of elements with abundances close to 
their NSE values at charged-particle reaction freeze-out. In 
particular, the entropies (~ 10 fcs baryon -1 ) and dynamical 
timescales (rd yn ~ 0.01 — 0.1 s) at freeze-out are similar to 
those explored in low entropy nucleosynthesis calculations 
by Woosley & Hoffman (1992). For example, in their model 
with S ~ 29 ks baryon -1 and Tdyn ~ 0.1 s, these authors 
find that all a— particles reassemble into heavy nuclei for 
Y e < 0.45, and that the final abundances closely resemble 
those of NSE (see the second column in their Table 2) . 

Neutron-rich NSE calculations were performed by Hart- 
mann et al. (1985). For Y e « 0.35 - 0.40, they find that the 



predominant isotopes created are 82 Se, 79 Br, and 78 Se (in 
order of decreasing Y e ; see their Fig. 3). The very low so- 
lar system abundances of these isotopes (total mass fraction 
Aq w 3 x 10~ 8 ) strongly constrains the rate of astrophys- 
ical events that produce material with Y e fa 0.35 — 0.4 in 
NSE: we estimate that the maximum allowed mass ejection 
of material with Y e ft! 0.35 - 0.4 into the ISM is M max » 
(M gal X /t gal ) « 1.5 x 1O" 7 M yr" 1 , where t gal w 10 10 yrs 
and M ga i w 5 x 10 10 Mq and are the age and baryonic mass 
of the Galaxy, respectively (Binney & Tremaine 1994). More 
neutron-rich material (Y e < 0.35) will produce even heav- 
ier isotopes and, possibly, an r-process. Since these elements 
are also very rare, a similar constraint can be placed on the 
rate that very low-Y e material is ejected into the ISM. For 
instance, isotopes near the second r-process peak (A ~ 130) 
have a typical solar mass abundance Xq ~ 10 -9 relative to 
hydrogen; therefore, if material ejected with 0.2 < Y e < 0.35 
produced equal numbers of N ~ 10 r-process isotopes in 
this mass range, the constraint on the mass ejected would 
be M max ~ 5 x 10~ 8 (A/10)Mq yr -1 , comparable to the 
rate for material with Y e pb 0.35 — 0.4 given above. In or- 
der to be conservative, below we use the constraint M max w 
1.5 x 10~ 7 Mq yr -1 for all material ejected with Y e < 0.4. 
However, a more precise determination of the isotopes pro- 
duced by the ejection of moderate entropy material with 
Y e ~ 0.2 — 0.35 would likely place even stronger constraints 
on the rate of CO mergers. 

Based on Figures 17191 and our one- zone calculations in 
32] we find that a fraction / a( j v ~ 0.2 — 0.5 of the disk's initial 
mass remains when it becomes advective at late times (see 
Table [TJ, with a higher fraction for less massive, more ex- 
tended disks. Of this material, we estimate that the majority 
of the mass will have Yj < 0.4, when averaging over many 
events with different initial disk properties (see Figs. 09] and 
Table [TJ. Although ejecta from mergers that occur outside 
Galactic disk will not enrich the ISM, population synthe- 
sis studies estimate that for large spirals such as the Milky 
Way only a small fraction (< 20%) of mergers occur in in- 
tergalactic space (Belczynski et al. 2006). Thus a fraction 
77 pb 0.1 — 0.5 of the disk's initial mass will both become un- 
bound from the disk and pollute the ISM with rare neutron- 
rich isotopes. Hence, if the average mass of the disks formed 
from CO mergers is (Md,o), their rate in the Milky Way 
cannot exceed 

Kr „ M max ln _ 5 f v N- 1 ( (Af d|0 ) V 1 _x 

where we have normalized r\ to a typical value ~ 0.2. 

Based on observed binary NS systems, Kalogera et 
al. (2004) find that the NS-NS merger rate in the Milky Way 
is between 1.7 x 10" 5 and 2.9 x 10" 4 yr -1 at 95% confidence. 
Population synthesis estimates (e.g., Belczynski et al. 2006) 
are consistent with this range but with larger uncertainties. 
Our constraint in equation (|17[l thus requires either NS-NS 
merger rates at the low end of current estimates, or a small 
average disk mass <C O.IMq. 

Ultimately, equation (I17p must be combined with GR 
simulations of the merger process that determine Md,o as a 
function of the total binary mass M to t = Mi + M2, the NS 
mass ratio q = M1/M2 ^ 1, and the NS equation of state 
(EOS). General relativistic merger simulations (Shibata et 
al. 2005; Shibata & Taniguchi 2006) find that when M tot 



© ???? RAS, MNRAS 000,rTHT2l 



Neutron-Rich Freeze- Out in Accretion Disks 11 



is above a critical value (j> 2.6Mq, depending on the EOS) 
the central object produced during the merger promptly col- 
lapses to a BH. In this case, Md,o decreases rapidly with 
increasing q (see Shibata & Taniguchi 2006, Fig. 13). For in- 
stance, when the double pulsar system PSR J0737-3039A/B 
(with q ~ 0.935 and M tot = 2.587M Q ) merges in w 100 
Myr, current simulations predict that the disk will have a 
mass ~ 10~ 3 — 10 _2 Mq if BH formation is prompt (Shibata 
& Taniguchi 2006). This disk mass is reasonably consistent 
with our constraint in equation (|17[) an d current NS-NS 
merger rate estimates. It is also consistent with the rela- 
tively low isotropic energies of short GRBs with measured 
redshifts (~ 10" 5 - 1O" 3 M c 2 ; Nakar 2007). 

On the other hand, simulations find that when M to t is 
below the threshold for prompt collapse, a hypermassive NS 
supported by differential rotation is initially formed (Baum- 
garte et al. 2000). If this NS is able to transport angular 
momentum outwards, it will collapse to a BH on a longer 
timescale (f» 100 ms; Shibata et al. 2006; Duez et al. 2006), 
producing a sizable disk (> 0.03Mq) whose mass is rela- 
tively independent of q (Oechslin & Janka 2006; Shibata et 
al. 2006). Since the disk masses formed after delayed col- 
lapse are comparable to those allowed by our constraint in 
equation (|17[) even for low TV ~ 10 -5 yr _1 , our results sug- 
gest that the formation of a long-lived hyper-massive NS is 
a rare outcome of NS-NS mergers. 

Since there are no known BH-NS binaries, the BH-NS 
merger rate is even less certain. Bethe & Brown (1998) argue 
that BH-NS mergers could be substantially more common 
than NS-NS mergers, with Bethe et al. (2007) estimating 
a rate ~ 10 4 Gpc _3 yr _1 , corresponding to ~ 10 -3 yr _1 in 
the Milky Way. Our results rule out such a high rate unless 
(Md,o) is less than ~ 10 _3 Mq. A low average disk mass 
may be possible, however, because only systems for which 
the NS-BH mass ratio falls within a fairly narrow range may 
produce a sizable disk upon merging (e.g., Miller 2005). 

The late-time nucleosynthesis in CO merger disks also 
places interesting constraints on the beaming fraction of 
short GRBs if most of these events truly result from CO 
mergers. Using the observed local short GRB rate of « 10 
Gpc 3 yr _1 (Nakar et al. 2006) and assuming that the merger 
rate is proportional to the blue stellar luminosity (Phinney 
1991), the short GRB rate in the Milky Way is estimated to 
be TVs grb » 10~ e yr" 1 (Nakar 2007). Thus, if all mergers 
that produce disks also produce short GRBs, the beaming 
fraction must obey 

This constraint implies that the jet half-opening angle must 
exceed ~ (2/ b ) 1/2 ~ 3O°(?7/O.2) 1/2 (A/ di o/O.lM ) 1/2 . This 
is consistent with other evidence that short GRBs may be 
less collimated than long GRBs (e.g., Grupe et al. 2006; 
Soderberg et al. 2006), which are instead inferred to have 
0.002 <: f b <: 0.01 (Frail et al. 2001; Bloom et al. 2003; 
Guetta, Piran, & Waxman 2005). Note that equation (I18|l is 
probably a stronger constraint than equation (|17|) because 
only CO mergers that form sizable accretion disks are likely 
to produce GRBs in the first place. 

In addition to material driven from the disk by late time 
winds, neutron-rich matter may be ejected dynamically dur- 
ing the merger process itself (e.g., Freiburghaus et al. 1999; 



Rosswog et al. 1999). Although the amount of material 
ejected in merger calculations depends sensitively on de- 
tails such as the NS spin and equation of state (Rosswog 
et al. 1999), the rare isotopes synthesized by this highly 
neutron-rich ejecta may provide a constraint on the CO 
merger rate that is comparable to that from late-time out- 
flows given by equation (117[) . 

Our conclusion that accretion disks formed from CO 
mergers freeze out neutron-rich is also consistent with the 
absence of SN-like longer wavelength transients coincident 
with these events. If CO merger disks froze out with Yj > 
0.5 instead of Yj < 0.5, their late-time outflows would 
efficiently synthesize 56 Ni (e.g., Seitenzahl et al. 2008), 
which could power an optical/infrared transient on ~ day 
timescales following short-duration GRBs. Despite inten- 
sive searches, however, such transients have not yet been 
observed. For instance, early optical follow-up of the GRB 
050509b (Hjorth et al. 2005) limited the amount of 56 Ni to 
< 1O~ 2 M (Kulkarni 2005; MPQ08), consistent with the 
modest amount of material with Y/ > 0.5 that we find in 
Figures 09] Although we do not expect much 56 Ni to be 
produced in late-time outflows from CO merger disks, ear- 
lier neutrino-driven outflows may produce up to ~ 10 _3 Mq 
in Ni, which could power a fainter optical or infrared tran- 
sient (MPQ08). In addition, outflows powered by fall-back 
accretion onto the BH at late times could power an X-ray 
transient on a timescale of days to weeks (Rossi & Begelman 
2008). 
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